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Abstract 

^ ■ Puiseux series are power series in which the exponents can be fractional and/or 

j — . negative rational numbers. Several computer algebra systems have one or more 

built-in or loadable functions for computing truncated Puiseux series. Some are 

,__, ' generalized to allow coefficients containing functions of the series variable that are 

dominated by any power of that variable, such as logarithms and nested logarithms 

C/3 ! of the series variable. Some computer algebra systems also have built-in or loadable 

c/3 functions that compute infinite Puiseux series. Unfortunately, there are some little- 

known pitfalls in computing Puiseux series. The most serious of these is expansions 
within branch cuts or at branch points that are incorrect for some directions in the 
£>. \ complex plane. For example with each series implementation accessible to you: 

^sO ' Compare the value of (z 2 + z 3 ) 3 ' 2 with that of its truncated series expansion 

about z = 0, approximated at z = —0.01. Does the series converge to a value that 

I/-) \ is the negative of the correct value? 

-^ ' Compare the value of ln(z 2 + z 3 ) with its truncated series expansion about 

^^ ■ z = 0, approximated at z = — 0.01 + 0. li. Does the series converge to a value that 

is incorrect by 2-7ri? 
• • ' Compare arctanh(— 2 + ln(z)z) with its truncated series expansion about z = 0, 

. £^ ■ approximated at z = —0.01. Does the series converge to a value that is incorrect 

^ , by about iri? 

At the time of this writing, most implementations that accommodate such series 
exhibit such errors. This article describes how to avoid these errors both for manual 
derivation of series and when implementing series packages. 
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1 Introduction 

This article is a companion to reference [5j. That article describes how to overcome 
design limitation that make many current Puiseux-series implementations unnecessarily 
inconvenient, such as not providing the order that the user requests or not allowing 
requests for negative or fractional orders. 

In contrast, this article describes how to overcome the more serious problem of results 
that are incorrect on branch cuts and at branch points for most current implementations. 
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This article is relevant to both truncated and infinite series of almost any type, includ- 
ing hierarchical, Fourier, Dirichlet and Poisson series. However, for concreteness the 
discussion is specific to Puiseux series that are generalized to permit in the coefficients 
sub-polynomial functions of the expansion variable, such as logarithms. 

Section [2] discusses branch bugs for logarithms. Section [3] discusses branch bugs for 
fractional powers. Section H] discusses branch bugs for inverse trigonometric and inverse 
hyperbolic functions. 

Let z be a complex variable, and let x, y, r, and 9 be real variables. An appropriate 
substitution can always transform any expansion point including oo, — oo and the complex 
circle at radius oo to z = x + iy = re td = 0. Therefore without loss of generality the 
discussion assumes that is the expansion point. 

To test an implementation for branch bugs in a truncated series U (z) for an expression 
u{z): 

1. Evaluate \u (0) — U (0)|. If there is a singularity at z = 0, then the result might be 
undefined even if the series is correct. Otherwise this absolute error should be or 
very nearly so. 

2. Do a high-resolution 3D plot of \u (z) — U {z) \ \ z — > x + iy centered at z = 
with enough terms to span an exponent range of at least 4. Try zooming in from 
a moderate initial box radius. An initial box radius of 0.5 works well for most 
examples in this article. 

3. 3D plots can easily miss discontinuities that are ribs, crevasses, or thin cusps ema- 
nating from z — - particularly if an edge doesn't lie along a grid line. Therefore 
if step 2 doesn't reveal an incorrect result, then 

(a) Do a 2-D plot of \u (z) — U (z)\ \ z — ¥ r e ie for 9 = (— it, n] and various fixed 
Tq that are well within the estimated radius of convergence. 

(b) For each critical direction 9 C defined in Section |2.1.1| plot \u (z) — U (z)\ \ z — > 
re l9c for r = [-R, R] and various fixed R > that are well within the estimated 
radius of convergence. 

4. If the result of step 1 is undefined because of a singularity, then exclude z = or 
clip the plot magnitude to a positive value ^ 1. 

5. Within rounding error and the radius of convergence, the surface and the curves 
should converge to 0.0 as the number of terms increases. If instead any of these 
plots converge to an obvious jump touching z = 0, then the formula is almost 
certainly incorrect. If there is a hint of a jump that grows from magnitude 0.0 at 
z = 0, then try instead plotting the relative error \{u (z) — U (z)) /u (z)\, excluding 
z = 0. 

For a real variable x, it suffices instead to evaluate -u(0) — U(0) and to plot \u (x) — U (x)\ 

and 

u(x) — U (x) 

u (x) 



2 Branch bugs for In 



"Spare the branch and spoil the child." 
- adapted from King Solomon's proverbs. 

Table [1] gives several logands together with one or more correct alternatives for their 
dominant 0-degree terms of the logarithm series expanded about complex z = or real 
x = 0. 

If an implementation gives a degree-0 term that isn't equivalent to these correct 
alternatives at x = and as x — > from both directions or at z = and as z — > 
from all directions, then the computer algebra result is incorrect. For example, most 
implementations currently give the following generalized infinite generalized Maclaurin 
series or a truncated version of it for In (z 2 + z 3 ): 

fc=0 

Within the radius of convergence 1, this series converges to values that are too large by 
2iri wherever 



y > A x < 
or too small by 2m wherever 

H<0Ai< 



1 + 3y 2 - 1 



1 + 3y 2 - 1 



with z = x + iy. For example, at z — O.li, ln(z 2 + z 3 ) ~ 4.6002— 3. 04192i, whereas 
series (JT]) truncated to o (z 4 ) gives approximately 4. 6002+3. 24126i. Series (JTJ is also 
incorrect at z = —0.1. 

If an implementation gives a degree-0 term that is equivalent to one of the results listed 
in Table [TJbut more complicated, then there is room for improvement of the simplification 
in ways described below. 

One way to avoid returning an incorrect result is to refuse attempting series expansions 
on branch cuts and on the branch points at their ends, either returning an error indication 
or an unsimplified result such as "series (• • • )". However, this precludes useful results for 
many examples of frequent interest, such as 

• fractional powers and logarithms of many expressions whose dominant exponent is 
non-zero or whose dominant coefficient isn't positive, 

• arcsin (u (z) ) , aiccos (u(z)), arccosh (u(z)), and arctanh (u(z)) at u(z) = 1 or 
u(z) = -1. 

Another way to avoid returning an incorrect result is to force the user to specify a numeric 
direction 8q for the series expansion variable z = re 10 , compute 



series 



ies (/ (re ie °) , r = + , o{r n )) 



substitute r — y ze~ l6 ° into the result, then preferably attach to the result the constraint 
"| z — re t0 ° A r > 0". The result is then guaranteed only for direction 9q. This is a 
reasonable approach when the only purpose of the series is to determine a uni-directional 
limit of an expression via that of its dominant term, and for a bi-directional limit we 
can invoke series (...) twice with two different values of 9q. However, this approach isn't 
appropriate for omni-directional limits. 

Moreover, with this approach it is important for # to have no default, such as the 
most likely choice 0. Otherwise many users won't realize that their formula might not 
be correct for non-positive or non-real z. 

In contrast, this article presents formulas that are correct for all 8q that aren't pre- 
cluded by any constraints provided by the user, such as "... | z > 0" or "... | — 7r/3 < 
argz < 27r/3". Moreover, even for the approach of requiring a numeric 8q, the formulas 
presented in the remainder of this article are helpful for determining the correct behavior 
for that #o- 

2.1 Incorrect extraction of the dominant term 

"The devil is in the details." 
- after Gustave Flaubert. 

Most computer-algebra systems use a particular branch when a multiply-branched 
function is simplified for numeric arguments. This branch is most often the principal 
branch. However, some computer algebra systems offer the option or the default of using 
the real branch for fractional powers having odd reduced denominators together with real 
radicands. Either way, for consistency the same branch should be used for expressions 
and their series. 

One source of incorrect In series is omitting the Ti term in the following universal 
principal-branch formula for the distribution of logarithms over products: 

\n(uv) = ln(u)+ln(i;)+Tz, (2) 

where 

T = axg{uv) — arg(w) — arg(t> ) (3) 

2tt if arg(w) + arg(t>) < — n, 

-2n if arg(-u) + arg(t> ) > n, (4) 

^ otherwise. 

This can be proved from 

In (|H) = hi(M)+ln(M), (5) 

ln(|iy|) = In (u>) — arg (w) i, (6) 

arg(ut>) = mods (arg(w) + arg(t>), 2n) . (7) 

Here mods {u^v) is the residue of u of mod v in the near-symmetric interval (—v/2, v/2] 
for v positive^ 



1 The name mods is inspired by that built-in Maple function. 



Remark 1 . These formulas require the useful but non-universal definition 

arg(O) := 0, (8) 

as is done in Mathematica® . If a built-in arg (. . .) function does anything else, then an 
implementer should prepend here and throughout this article appropriate cases for each 
possible combination of an argument of arg (. . .) being 0. For example, 



T 



if u = V v = 0, 

arg(m>) — arg(u) — arg(t>) otherwise. 

ifu = 0Vf = 0V — tt < arg(u) + arg(w) < tt, 

2tc if arg(w) + arg(u) < — tt, 
-2tt otherwise. 

Here and throughout this article, Boolean expressions and braced case expressions are 
assumed to be done using short-circuit evaluation from left- to-right within top-to-bottom 
order to avoid evaluating ill-defined sub-expressions and to avoid the clutter of making 
the tests mutually exclusive. 

Alternative (|3J) is more compact than alternative (J3J) and reveals that jumps in T can 
occur only where one of arg(w), arg(t> ) or arg (uv) is tt. However, alternative (jlj) is more 
candid because it makes the piecewise constancy manifest rather than cryptic. Moreover, 
approximate values are often substituted into expressions for purposes such as plotting, 
and the conditional alternative (jl]) avoids having the magnitude of the imaginary part of 
a result be several machine e when it should be 0: Unlike the unconditional alternative, 
the conditional alternative never subtracts two approximate angles from approximately 
tt, giving approximately 0. 

For series (In (...), z — 0, o (z n )) with negative n, the result is + o (z n ) if the logand 
doesn't contain an essential singularity. In contrast, for a non-negative requested n, the 
usual algorithm for computing the logarithm of a series entails converting the dominant 
term of the logand series to 1 by factoring out the dominant term then distributing the 
logarithm over the resulting product: 

In (c (z) z a + g(z)) -> In (c (z) z a (l ' <A ' '' 



c{z) z c 



-► (^ + m( C (^))+lnhH--^Lj. (9) 

Here c (z) z a is the dominant term and g(z) is the sum of all the other terms, with 



Many of the formulas involving arg (...) in this article can be expressed more concisely using the 
unwinding number described in[lj and later redefined more conveniently with the opposite sign in [2]. 
But alas, unwinding numbers aren't yet built-into the mathematics curriculum and most computer- 
algebra systems. 



alternatives ([3]) and (j3j) giving 

Q = arg (c(z)z a + g(z)) - arg (l + 4^) ~ arg (c (z) z Q ) (10) 

2vr if arg M + ^^J + arg (c (z) ^ Q ) < — k, 

-2tt if arg (l + ^L J + arg (c (2) ^ Q ) > vr, (11) 

otherwise. 

Always Q = at z = 0. 

It is important to simplify Q as much as is practical for each particular logand series. 
Not only is the result more intelligible - it is usually also more accurate for approximate 
computation: There can be catastrophic cancellation between terms of g(z). Therefore 
without good algebraic simplification, Q can be dramatically incorrect along and near 
branch cuts when evaluated with approximate arithmetic. 

Remark 2. For example, Q = if g(z) = 0, because then arg(l + g(z)/ (c(z)z a )) = 
and arg (c (z) z a ) must be in the interval (— 71", n]. 

Even if we can't determine a simpler formula that is equivalent to formula (ITT]) 
throughout the entire complex plane or the entire real line for real z, Q might be equiva- 
lent to a simpler expression u throughout the radius of convergence or the useful portion 
thereof. It is especially important to exploit this for examples such as 

series (z^ 1 + e z In (-2 - z) , z = 0, o (z 5 )) 

where O or u infects all but one result term, giving a bulky result that is difficult to 
comprehend. 

Let R > be the classic radius of convergence computed disregarding any closer 
branch cuts, or let R be the "radius of computational utility" for divergent series. Let 
R > be the largest radius from z = within which Q and a simpler u give identical 
values. Radius R can be an arbitrarily small portion of R because a branch cut can pass 
arbitrarily close to z = 0. However, it is almost always justifiable to use u in place of Q 
because: 

• If our purpose is to determine the local behavior of the series at z — 0, such as for 
computing a limit by computing the limit of the dominant term, then any R > 
is sufficient justification for using u. 

• It seems pointless to use Q rather than a simpler u for any purpose if R > R or if 
R is greater than the percentage of R beyond which convergence is impractically 
slow or subject to unacceptable catastrophic cancellation. 

• If we don't expect a generalized Puiseux series to capture infinite magnitudes as- 
sociated with singularities not at z — 0, then why should we expect such series 
to capture the less severe finite-magnitude jumps associated with branch cuts that 
don't touch z = 0? 

• We can take the view that the generalized radius of convergence is the distance to 
the nearest singularity or jump in Q that we can't account for with u, and there 
should be no expectation that a series is truthful beyond its generalized radius of 
convergence. 
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Here is one such opportunity for computing an u that is significantly simpler than Q: 

Proposition 3. Let u> rca \ denote u> for the special case of real z. If all of the terms in 
the truncated logand series have real coefficients and integer exponents, then u Tca \ = 0. 

Proof. If the truncated logand series has all integer powers of real z and all real coeffi- 
cients, then 1 + g (z) / (c (z) z a ) is real for all real z, and arg (1 + g (z) / (c (z) z a )) = for 
all real z such that g (z) / (c (z) z a ) > — 1. There is a singularity wherever g(z)/ (c(z)z a ) = 
— 1, providing an upper bound on the radius of convergence. Also, — it < arg (c(z)z a ) < n. 
Therefore throughout the radius of convergence 

-7T < arg M + ^j4^j + arg (c (z) z a ) < n 

in equation (fTTj) . making w rca i = 0. D 

Remark 4. If we are in a mode that consistently uses the real branch for fractional 
powers having odd denominators, such as the TI-Nspirdj real mode, then more generally 
cereal = if all of the coefficients are real and none of the reduced exponents have even 
denominators. 

2.1.1 The non-zero dominant exponent case 

"Beware of geeks bearing formulas." 
- Warren Buffett. 

Here is a useful easy simplification test for real or complex z when the dominant 
exponent q^O: 

Proposition 5. If all of the coefficients in a truncated logand series U are real and all 
of the exponents of z in U are integer multiples of a non-zero dominant exponent, then 
w = 0. 

Proof. Let the truncated logand series be U = c (z) z a + g(z) with dominant term c (z) z a 
having non-zero a. Also, let "near z = 0" denote all 



re %e I r>0 A 



g{z) 



c(z)z a 



<1, (12) 



which bounds the radius of convergence. From equation (llOl) . it is apparent that fl can 
jump only where one of its three arg (. . .) terms jumps. Term arg (1 + g(z)/ (c(z)z a )) 
can't jump near z = because c(z)z a dominates g(z). Every term in c(z)z a + g(z) 
can be expressed as k(z) (z a ) m with real k(z) and integer m, because the transformation 



2 The computer algebra used in Texas Instruments products has no name separate from the variously 
named calculators, Widows and Macintosh products that contain it. The most recent such product is 
named TI-Nspire. 



z 1 — > (z a ) m with ?7i = 7/a is always valid for integer m, as discussed in Section 13.11 



Thus c (z) z a + g(z) is real where 



arg (c (2; J 2 



7T, 



(13) 



which are the only places that arg (c (z) z a ) jumps. Since c (z) z a is negative there and 
dominates g(z), adding the relatively small-magnitude real g(z) to negative c (z) z a leaves 



axg(c(z)z a + g(z)) = n 



(14) 



near z = 0. In equation (llOp . arg (c (z) 2:° + g{z)) and arg (c (z) z a ) have opposite signs. 
Therefore these jumps cancel within the radius of convergence. □ 

Propositions [3] and O don't identify all opportunities for dramatically simplifying oj. 
Consider equation (jlOp in the neighborhood of z = 0. Let 

T)(8) = lim arg (c (re l6 ) ) . (15) 

r->0+ 

Then in the punctured neighborhood of z = 0, the term arg (c (z) z a ) is 

lim arg (c(re ie ) (re ie ) a ) = mods {n (9) + a9, 2vr) . (16) 

i — >o+ 

Let fj denote 77 (0) in formula (TTB]) in the common case when 77 (#) is independent of 
9 for all — 7r < 9 < -k not precluded by constraints provided by the user or introduced by 
the computer algebra system. Solving mods (?) + a9 c , 2tt) = ir for the critical angles 9 C 
gives 

(On -I- 1 W _ f> 

(17) 



(2n + l)7T-77 
7T < r = ± —^ < IT, 



a 



with integer n. Solving the inequalities in (TTTj) for n gives for a > 



versus for a < 



77 

_2tt " 


1 a 

~ 2 ~ 2_ 


< n < 


77 

_2tt " 


1 a 
" 2 + 2_ 


2ir ' 


1 a 
-2 + 2 


< n < 


2tt 


1 a 

2 2 



(19) 



There are no such angles if |a| < 1 and 7} is sufficiently close to 0, in which case 
00 = 0. More generally let 77 be a lower bound and fj be an upper bound on 77 (9) over 
— 71 < 9 < 7r not excluded by any constraints. Then there are no critical angles if 

(|a| -1)tt <?7 A 77 < (1- |a|)7r. (20) 

Here is how we can proceed when none of the above tests are beneficial: As z — > 0, 

9(z) 



c(z)z a + g(z) = c(z)z a [ 1 + 



c(z)z a 



— ¥ c[z)z 



(21) 



Therefore the solution curves to equations (fT3|) and ( |T4"|) pair to form cusps where 
they don't coincide to cancel. Along a critical angle 9 C , arg (c(z)z a ) = it. If also 
arg (1 + g{z)/ (c(z)z a )) = along 9 C near z = 0, then 9 = 9 C is also the companion 
solution to equation (|14|) near 2 = 0. The arg (. . .) terms for equations (fT3|) and (J14J1 
have opposite signs in formula f lTOj) . so these jumps cancel within the radius of conver- 
gence. Since arg (1 + g(z)/ (c(z)z a )) can't jump near z = 0, u is then identically in the 
angular neighborhood of 9 C near z — 0. 

If u; = in the angular neighborhood of every critical angle, then u = 0. For 
example, this is true for In (— z 2 + r 3 ) and In (— z -2 + z), which don't satisfy Proposition 
[5] or inequality ([20]) . Even if there are some non-zero cusps, constraints might exclude 
those cusps for a non-infinitesimal distance from z = 0. Even if there are some included 
cusps, it might be possible to omit either the 2ir case or the — 2n case as follows: 

For a > 0, as 9 increases through 9 C , arg (c ■ (re 10 ) ) increases with 9 on both sides 
of a jump down from it to (— n) + . If arg (1 + g(z)/ {cz a )) > along 6 C near z — 0, then 
from equation (fTT|) . u = —2n along and clockwise of 9 C until but excluding the curved 
solution to equation (I14p . If instead arg (1 + g(z)/ (cz a )) < along 9 C near z — 0, then 
to = 2ir counter-clockwise of 8 C through the curved solution to equation (|14j) . 

Similarly for a < 0, for which arg (c • (re l9 ) ) decreases on both sides of a jump up 
from (— 7r) + to 7r as 9 increases through 9 C : If arg (1 + g(z)/ {cz a )) > along 9 C near 
z = 0, then u; = — 2tt along and counter-clockwise of 9 C until but excluding the curved 
solution to equation ( JT4"|) . If instead arg (1 + g(z)/ (cz a )) < along 9 C near z = 0, then 
u) = 2ic clockwise of 9 C through the curved solution to equation (fl4l) . 

Thus we can omit the 2n case if for all critical directions arg (1 + g(z)/ (cz a )) > 0, 
or we can omit the — 2ir case if for all critical directions arg (1 + g(z)/ (cz a )) < 0. For 
example, we can omit the 2tt case for In (iz 2 + iz A ) and we can omit the — 2tt case for 
In (z + iz 2 ). 

Notice that although u can be 2n either clockwise or counter-clockwise of a critical 
direction, u can't be 2ir along a critical direction. 

For real t let 

1, if t < 0, 
signum(t) := { 0, if t = 0, 

if t > 0. 

To determine whether expression arg (1 + g(z)/ (cz a )) is zero, positive or negative along 
9 C near (but not at) z = 0, we could attempt computing 

g ire 19 ") 
lim signum arg 1 + 




>o+ \ \ c • [re 1 ' 

However, such limits are beyond the capabilities of most computer algebra systems if g 
has more than one term - particularly if any of the exponents are fractional and/or any 
of the coefficients depend on z. A simpler surrogate within the radius of convergence is 
to instead compute 

lim signum (^(^4)). (22) 



However, this limit is also often beyond the capabilities of most computer algebra systems. 
Fortunately there is an easily-computed surrogate for formula (I22p : Let s (re* 9 ) be the 
coefficient and (3 be the exponent of the lowest-degree term of g(z)/ (c(z)z a ) for which 

d( z ) n\._ i- ex f„f„J6. 



cz c 



6A := lim + 3(s(re iec )e i/3ec )) + 0, (23) 



if any such term exists. Let I c (. . . , 6 C ) := if no such s (re' 61 ) exists. The imaginary part 

of s (re l9c ) (re l6lc ) dominates the imaginary parts of all subsequent terms of g(z) along 
9 C if any exist. Therefore I c is a more-easily computed surrogate for determining whether 
arg(l + g(z)/ {cz a )) is zero, positive or negative along 9 C near z — 0. The term limits 
for computing I c are trivial in the common case when a coefficient is a numeric constant, 
and easy even if a coefficient is a typical sub-polynomial function of z. We can use Q 
given by formula (fTUj) or (TTTj) if a coefficient contains an indeterminant other than z that 
isn't sufficiently constrained to decide the sign of the imaginary part and if no imaginary 
parts were non-zero for previous terms. 

The order to which the logand series is computed might not reveal the lowest-degree 
term for which one of the I c (. . . , 6 C ) isn't 0, thus affecting the resulting value of Q. For 
example, if the series for In (z 2 — iz 3 — iz 6 ) is computed to o (z n ) with n > 4, then the 
dominant term is 

, / 2 \ J ^ n ar S (I — iz — iz 4 + ■ ■ ■ + o (z n )) + arg (z 2 ) < — n, 

111 \Z j \ \ 

I otherwise. 

However, the piecewise term is absent if the series is computed to o(z 3 ), making it 
incorrect by about 27ri for values such as z = z/10 — 10 -6 and z = — z/10 + 10~ 6 . It is 
disturbing that computing additional terms of a logand can thus change the zero-degree 
term of the logarithm series. Moreover, it is awkward for algorithms that compute 
additional terms incrementally, such as described by Norman [4J. Such is the nature of 
partial information along or near a branch cut. 

A way to avoid this annoyance is to use the original logand expression in the arg (...) 
sub-expressions of alternative ( TTOj) or ( TTTj) rather than using a truncated series for that 
logand. However, the original logand isn't always available: Perhaps as the logand we 
are given a truncated series, not knowing a closed-form expression that it approximates. 
Even if we did know, using the original expression can make computing a series in one 
step give a different result than series composition. Such compositional inconsistency is 
an undesirable property. 

More seriously, including the original non-series expression as a proper sub-expression 
of a series result is worse than simply returning the original expression rather than a series: 
Presumably the user requested the series to obtain insight about the behavior of the 
function, or to enable symbolic operations that otherwise couldn't be done, or to enable 
a numeric approximation. Returning a result that contains the original expression as a 
proper sub-expression thwarts all of these objectives. The annoyance of order- dependent 
coefficients is vastly preferable. A consolation is that a zero-degree term that differs so 
dramatically from that of the infinite series is at least appropriate for a nearby problem 
across the branch cut. 
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Remark 6. If the truncated logand series is simple enough so that for remaining cusps 
we can solve 

arg (c (x + iy) (x + iy) a + g (x + iy)) = n (24) 

for x (y) or for y (x), or else solve 

arg (c (re ie ) (re ie ) a + g (re 10 )) = n (25) 

for 9 (r) or for r (8), then we can construct a more candidly explicit representation of oj. 
For example, with In (z 2 + z 3 ), 



n 



2tt if3(z)<0A0<SR(z)< 

-2tt if 3 (z) > A < 9? (2) < 
otherwise. 



y/l + 3^) 2 - 


-1 


3 


1 


\A + 3^) 2 - 


- 1 



(26) 



There is catastrophic cancellation here computing W 1 + 33 (2) — 1 for |£$ (z)\ C 1. 

This can be avoided by expanding the expression into the series 3 (z) /2 — 33 (z) /8 + 
. . ., which is also more consistent with the user's request for a series result. Perhaps 
generalized series reversion could be used to obtain an explicit truncated series solution 
to equation ( 12~4|) or (j2"5"|) even when we can't solve them exactly - at least when the 
coefficients are all numeric and the exponents are all non-negative integers. 

Remark 7. When z is real, we can almost always simplify u rea \ to either or a two-piece 
result that is — 2n for negative z and/or positive z, but everywhere else as follows: 
If isn't a critical angle or I c (g(z)/ (cz a ) , 0) < 0, then u; rca i = for z > 0. Otherwise 

Wreal = — 27T for Z > 0. 

If 7r isn't a critical angle or I c (g(z)/ (cz a ) , vr) < 0, then u rc& \ = for z < 0. Otherwise 

Wreal = ~ 2vT for Z < 0. 

2.1.2 The degree-0 case 

If a = and r/ < 71", then u; = because for formula (II ip . arg (1 + g'(^)/ (c(^)z a )) — > 0. 
If instead a = and arg (c (z)) = 7r throughout all z near 2; = that aren't excluded by 
any constraints, then arg (1 + g{z)j (c(z)z a )) + arg (c (z)) can't be less than — ir, because 
arg(l + g(z)/ (c(z)z a )) — > 0. Therefore we can at least omit the 2tt case. Thus using 
00^^ to denote u for this special case of z = being in the branch cut, 

[-2* if arg f 1 + 44) > 0, r97 , 

w ff)0 = < V. c(z)J (27) 

^0 otherwise. 

There are singularities wherever g(z s )/c(z s ) = —1, and those singularities aren't at 
the expansion point z = 0. Consequently the radius of convergence doesn't exceed the 
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least of those \z„\. Therefore 

f— 2tt ifO<arg( 44 I <*, 
^,o = { \c(z)J (28) 

otherwise. 

Formula f l28l avoids adding a small magnitude quantity to 1.0 near z = 0, so this formula 
is more likely to be correct than formula (1271) with approximate arithmetic. 

Whenever we are comparing arg (/ (z)) with — 7r/2, 0, 7r/2, or 71", we can often simplify 
the test further and make it more accurate for approximate numbers by writing the 
comparison in terms of dt(f (z)) and/or Q (f (z)) - particularly when / (z) has more 
than one term. Thus 

oJn,o { \c(z)J (29) 

otherwise. 

This test can often be simplified further: Let b(z)z cr be the dominant term oi g(z)/c(z), 
and let 

r(0) = lim arg (b (re ie )) . (30) 

r— >0+ 

Let f denote r in the common case where it is independent of 9 for all < 9 < n that 
aren't excluded by any constraints. 

In the neighborhood of z = 0, as 9 increases from (—it) through tt, arg (g(z)/ (c(z)z a )), 
arg (b(z)z cr ) and arg (fz a ) increase from (-<T7r) + + f through air + f , but jumping down 
from 7r to (—7r) as 9 increases past every critical angle where mods(a9 c + f, 2ir) = it. 
Expression Q (g (z) / c(z)) changes sign at those places and also at every critical angle 
where mods(<7# c + f, 2tt) = 0. There are no critical angles of either type if < a < x ji 
and r is sufficiently close to 7r/2 or —it/2. More specifically if r is a lower bound on r(9) 
and r is an upper bound, then co^o = when 

(a-l)n<TAr< -air. (31) 

If instead 

an < t A T < (1 - a) 7T, (32) 

then 

^,o = < . (33) 

I — 2tt otherwise. 

As examples, w ff) o = for In (—1 + 22; 1//4 ) , whereas equation fl33|) applies to In (—1 — zz 1 / 4 ) . 

Remark 8. If critical angles occur only at the edges of regions not excluded by constraints, 
then we can compute I c (g(z)/c(z), . . .) at those critical angles and at one included non- 
critical angle to attempt simplifying formula (129|) . For example with In (—1 + iz 1 ^ 2 + z), 
the one critical angle is n, along which I s is non-positive. I c is also non-positive along 
the included non-critical angle —tt/2, so oj nt0 = 0. 
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In contrast for In (—1 — iz 1 ! 2 + z) , the one critical angle is ir, along which I c is non- 
positive. However, I c is positive along the included non-critical angle n/2. Therefore u 
is for the ray z < but — 2tt everywhere else. 

If a critical angle is interior to an included region, then we can simply use formula 
fl29|) . Remark [7] is also applicable to uj^a for real z. 

Remark 9. For an example such as In (c + z 2 + z 3 ) where a value such as 1, -1 or could 
subsequently be substituted for the literal constant c, a completely correct result should 
piecewise account for all possible cases. 

2.2 Incorrect extraction of the dominant coefficient 

A logarithm of a product is more concise and more efficient to approximate numerically 
than a sum of logarithms of the factors. For this reason, users often prefer a result that 
contains a logarithm of a product rather than an equivalent sum of logarithms. 

However, at this time many existing generalized Puiseux-series implementations al- 
ways distribute the logarithm of the dominant term over its coefficient and cofactor. 
This distribution is justified for a hierarchical series if the coefficient is itself a series in 
logarithms or nested logarithms depending on z - at least when this distribution uni- 
fies two otherwise different logarithms. For example, it is more consistent to have In (z) 
throughout a series than to have In (z) in some places and In (— 2z) in other places. 

However, this distribution provides an additional opportunity for incorrect results 
caused by ignoring the Ti term in equation (T5]). Applying this transformation twice and 
simplifying gives 



\n(c(z)z a + g(z)) 



-» 



($i + In (c (z)) + In (z a )) + In 1 + 



9W 



c (z) z c 



(34) 



where 



$ = fi + arg(c(z)z a ) — arg(c(z)) — &rg(z a ) , 

{2n if arg (c(z)) + arg (z a ) < —n, 
—2tt if arg (c(z)) + arg (z a ) > tt, 
otherwise, 



arg (c (z) z a + g{z)) - arg I 1 + 



9W 



c (z) z' 
2tc if arg (c(z)) + arg (z a ) + arg I 1 + 

— 2tt if arg (c(z)) + arg (z a ) + arg I 1-f 



arg(c(z)) -argO a ) 
9{z) 



c{z)z' 
9{z) 



c[z)z c 



<-7T, 
>7T, 



(35) 
(36) 

(37) 
(38) 







otherwise. 



Alternatives (1351) and (1361) are advantageous when u is a constant, because only one 
term occurs in the arguments of arg (...). Otherwise alternative (1381) is more concise and 
candid. 
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Let v be a lower bound and v be an upper bound on mods (a8, 2tt) for all angles 
—tt < 6 < tt that aren't excluded by constraints, with rj and r\ being corresponding 
bounds on r\ (8) defined by equation (j!5p . In equation (155]) the 2tt case can be omitted if 

V + v > -7T, (39) 

and/or the — 2tt case can be omitted if 

rj + v < ix. (40) 

The same is true for alternative (I38p because arg (1 + g(z)/ (c(z)z a )) — > 0. 

We can use r\ — fj — fj for the common case where rj is independent of 6. For 
unrestricted complex z we can use v = min (1, |a|) vr and v = —v. For unrestricted real 
z we can use the intervals 

' [an, 0] if - 1 < a < 0, 

r -i I [0,0] if« = 0, , A1 . 

\v,v] = < (41) 

L J I [0,ott] if0<a<l, 

k [— 7r, 7r] otherwise. 



Alternative (1371) reveals that jumps in u near z = can occur only where arg (z a ) = tt, 
where arg (c(z)) = tt, and where arg (c (z) ^ a + 5 , (^)) = tt. When arg {c(z)) ^ 0, the three 
solution sets for these three equations typically don't pair to cancel or form cusps. Instead 
they typically form non-cusp wedges of values — 2tt, and 2tt. Therefore, if there are 
critical angles interior to the included directions, then: 

• Non-zero local (f) occur in infinitely more directions than non-zero u. Therefore it 
is less likely that any constraints will preclude all of the non-zero wedges or even 
all of the positive ones or all of the negative ones. 

• It is impossible for both edges of every wedge to coincide and thereby cancel to 
give = for all directions near z = 0. 

Thus we often pay dearly for distributing the logarithm over a non-positive coefficient 
and its co-factor, which is unnecessary in many applications. 

Remark 10. For real z, lim 2 _j. + ar g (-2°) — and lim z _> ar g (1 + g{z)/ (c(z)z a )) — > 0. 
Therefore equation (138]) reveals that = for z > if \im z ^.Q+ axg (c(z)) ^ {— tt,tt}. 
When lim z ^, Q + a.ig(c(z)) G {— tt,tt}, then we can compute I c (g(z)/ (c(z)z a ) , 0) from 
definition (1231) to decide whether <ft is 2tt, —2tt or for z > 0. Similarly lim 2 ^ - ar g (^ a ) = 
mods (air, 2tt). Therefore if 

lim (mods (air, 2tt) + arg (c(z))) fi {—tt, tt} , 

z— >0~ 

then = for z < 0. Otherwise we can compute I c (g(z)/ (c(z)z a ) ,tt) to determine 
whether <fi is 2tt, —2tt or for z < 0. We can then combine these results to obtain rea i- 
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2.3 Incorrect extraction of the dominant exponent 

If the dominant exponent is neither nor 1 and either the coefficient of the dominant 
term is 1 or we have distributed the logarithm over it, then we can consider extracting 
the exponent of the dominant power. 

2.3.1 Principal branch 

The relevant universally correct principal-branch identity for extracting the exponent of 
In (z a ) is 

\n(u a ) = aln(u) + fi, (42) 

where 

£ = arg(u Q ) -aarg(w), (43) 

= 2tt l- Qarg(M) . (44) 

[2 2tt J K J 

These alternatives can be derived from identity (^ together with ln(|-u a |) = aln(|w|) 
and the fact that arg (u a ) = mods (a arg (u) , 2ir). Alternative (144")) is more accurate 
for approximate arithmetic and more candid because it is manifestly piecewise-constant 
integer multiples of 2tt. 

Expression £ can be simplified to a single constant 2mr if 

(2n-l)n < ctarg(w) < (2n+l)7T, (45) 

for an integer n together with the smallest and largest arg(w) in (—71,71] that aren't 
excluded by any constraint. For example, £ = if 

— 71 < a arg(w) < it. (46) 

If instead (2n — 1)tt < ctarg(-u) < [2n + 3)tt for all included arg (u), then £ can be more 
candidly represented as 

2mr if ctarg(u) < (2n + 1)tt, 

[2n + 2)tt otherwise. 

Extraction of a fractional power has the benefit of converting a troublesome fractional 
power into a benign multiplication by a fraction. 

Also, there is more justification for extraction of the dominant exponent for a hierar- 
chical series wherein the coefficient can be a series in a logarithm or nested logarithm of 
z. 

In our case, u — z, and combining this transformation with distribution over the 
dominant term and the coefficient thereof gives the overall transformation 

\n(c(z)z a + g(z)) -> (m + In (c (z)) + a In (z)) + In ( 1 + -^-\ (48) 

V c(z)z a ) 
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where 



* 



$ + arg (z a ) — a arg (z) , 

1 a arg(z) 

2 2tt 



$ + 2tt 



f2 + arg (c(z)z a ) — arg (c(z)) — a arg (z) , 

27T if arg (c(z)) + arg (z a ) < —n, 



Q + 2-k 



1 a arg(z) 
2 



2tt 



arg (c (z) z Q + g(z)) - arg 1 + 



+ ^ — 2n if arg (c(z)) + arg (z a ) > it, 
otherwise, 



c (z) z c 



arg (c (2)) — a arg (z) 



2^ 



1 a arg(z) 
2 



2^ 



27r if arg(c(z)) + arg(z a ) + argl 1 + 



5l^J 



— 27r if arg(c(z))+arg(z a ) + arg( 1-f 
otherwise. 



c{z)z" 



c{z)z° 



<-7T, 



>7T, 



(49) 

(50) 
(51) 

(52) 
(53) 

(54) 



Let if) denote a perhaps simpler local version of \1/ near z = 0. Alternative (ISDj) is 
advantageous when ^ is non-constant but is simple. Alternative (|5"2"|) is advantageous 
when ^ and are non-constant but uj is simple. Otherwise alternative (|54|) is most 
candid and least subject to catastrophic cancellation. 

Equation (153|) reveals that the critical angles are the union of the solutions to arg (z) = 
7T, a,Tg(c(z)z a ) = it, and (if c(z) depends on z) arg(c(z)) = n. The piecewise-constant 
values of \l/ are also multiples of 27r, but no longer limited to be one of — 27r, or 27r. 
Therefore with or without constraints on arg (z) , the local version if) is even less likely 
than to be identically or otherwise constant, and likewise for if> iea \- However, we 
will see that the use of \1/ rather than Q is mandatory for the fractional power of a 
series. Therefore it is important to simplify \1/ as much as is practical: The floor term 
in alternatives (I50p . (I52p and ( 154"|) can be simplified as discussed at the beginning of 
this sub-subsection. Subsection 12.21 describes how to simplify the conditional term in 
alternative (152|) . Simplification of the conditional term in alternative (154]) is similar 
because arg (1 + g(z)/ (c(z)z a )) — > 0. 

For real z, we can compute a local version of \1/, if) Tea \, by computing rea i, then adding 



2tt 



a 
2 



to the case for negative z. 



2.3.2 The real branch for In of fractional powers 

The Derive computer algebra system offers a Branch control variable that, if assigned the 
value Real causes fractional powers having odd denominators to use the real rather than 
principal branch when a radicand is negative. TI-Nspire bundles this choice into its Real 
mode. 
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For example, in Real mode (—1) — > — 1 rather than 1/2 + i\/3/2, and (—1) — > 1 
rather than —1/2 + iy/3/2. This default option is much appreciated by students and 
faculty who fear or loathe non-real numbers, which is a majority of TI's customers most 
of the time. However, this choice is incompatible with formulas ( l4"2"|) through (I43p . For 
example, these principal-branch formulas give 



In (x 2/3 ) -> - In (x) + ( arg (x 2/3 ) - - arg (x) J 



i. 



If we subsequently substitute -8 for x, then this result simplifies to 

^ln(-8) + ^arg((-8) 2/3 )-^arg(-8)')^ -> 2ln(2) + -ttz, 

whereas for the real branch, In ((— 8) 2 ' 3 ) — > In (4) — > 2 In (2) . 

The easiest way to simplify logarithms of fractional powers using the real branch is to 
refrain from extracting reduced fractional exponents having an odd denominator unless 
the sign of the base can be determined - perhaps with the assistance of a constraint on 
arg (z) . Moreover, for real-branch mode we should also refrain from extracting integer 
exponents that are even. For example, we shouldn't extract the exponent 12 from In (x 12 ) 
for a real x because 

ln(x 12 ) -> 12 In (a;) + (arg (x 12 ) - 12 arg (x)) i 

, . , f Yl-ni if x < 0, 
-)- 12 In {x)-{ 

I otherwise. 

Even though this result is real for all real x, this appearance of i in the result is un- 
welcome to most real-mode customers. Even principal-branch complex-mode customers 
would rather not see i in an expression if it can as concisely or more concisely be ex- 
pressed without i - especially if the result is real for all real values of interest for any 
variables therein. However, to help reduce the number of distinct logands in a result for 
x G M, we can extract all of an even exponent but 2: 

ln(x 12 ) -)■ 61n(x 2 ). 

Another alternative for real x is ln(x 12 ) — > 121n(|x|). However, absolute- value func- 
tions are troublesome and best avoided where possible. For example, if in the same series 
we use 

, /ox , , , I 2iri if x < 0, 

In (x 3 ) ->■ 3 In (x) - < 

I otherwise, 

then we obtain a series that contains both In (x) and In (\x\). 
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3 Branch bugs for fractional powers 

Some series implementations can give incorrect series results for fractional powers. For 
example, with (z 2 + z 3 ) most implementations currently incorrectly give the equivalent 

V 3(-l) fc (2fc)! fc+3 

f^ Q (2k - 3)(2k - l)k\H k 

or a truncated version of it. This series is incorrect by a factor of -1 left of x = 
( \/3y 2 + 1 — 1 j /3 for z = x + iy. 

Table [2] lists one or more correct alternative multiplicative correction factors for the 
Puiseux series of some fractional-power expressions expanded about complex z = or 
real x = 0. 

The probable causes of an incorrect result are related to those for logarithms of 
series. The usual algorithm for computing a numeric power /3 of a series requires that 
the dominant term be 1. Therefore if the series of the radicand is c(z) z a + g(z), with 
c (z) z a being the dominant term, we must: 

1. Factor out c (z) z a . 

2. Distribute the exponent (3 over the three factors (z) , cz a and 1 + g (z) / (c (z) z a ). 

3. Transform (z a ) to z a ^. 

4. Compute the series for (1 + g (z) / (c (z) z a )) by the usual algorithm. 

5. Distribute c (z) z a/3 over the terms computed in step 5. 

Steps 2 and 3 can contribute to an angular rotation factor of the form (— 1)*" that should 
also be distributed in step 5. 

3.1 Principal-branch series of fractional powers 

Zippel's formula [6J generalizes to the following universal principal-branch formula for 
the distribution of real exponents over products: 

{uvf = (-lfuV (55) 

where 

5 = — (arg f (uv) ) — arg(-u /3 ) — arg(f /3 )j . (56) 

For real a and /3, a universal principal-branch formula for transforming a power of a 
power to an unnested power is 

(w a f -> (-lfw a P (57) 

where 

C := — (arg (w a ) — a arg (w)) . (58) 

7T 
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This can be derived from the identities 



|p| = (-i)- ar s(p)/> 

I a 1/8 i i a/8 

Applying these formulas to our case gives 

(c(z)z a + g(z)'f -> (c(z)z a (l + 



(59) 
(60) 



g{z) 

c (z) z c 



-)• (-l) A c(z/^(l + 



g(*) 

c (z) Z a 



(61) 



where modulo 2, which suffices for (— 1) A , 

A = i( &Tg (c(z)z a + g(z))-aj:g(l + ^- c 

71" \ \ C{Z)Z C 

IT 

with ty defined by alternatives (|49p through (}54"1) . Let 

L := (-1) A . 



arg(c (z)) — a arg(z) 



L1L2 



(62) 



(63) 
(64) 



where 



2/3 



L 2 



1 a arg (z) 
2~ 



2tt 



■1) 

'— 1) 2/3 if arg(c(z)) + &rg(z a ) < —n, 

I i '-' '"^ <j {-l)- w if arg(c (z)) + arg(z a ) > rr, 

1 otherwise, 



(65) 
(66) 



•1) 2/3 if arg(c(z)) + arg(z") + arg I 1 + 



g{z) 



C[Z)Z C 



■1) 2/3 if arg(c (z)) + arg(z Q ) + arg I 1 + 



9W 



c(z)z c 



< -7T, 



> 7T, 



(67) 



1 



otherwise. 



Alternative (166]) is preferable for L 2 when u is constant. Techniques described in 
subsection 12.31 can further simplify the floor sub-expression in L\. Also, the exponents 
of -1 in these formulas can be simplified and canonicalized by replacing them with their 
near-symmetric residue modulo 2. 

Remark 11. In the common case where (3 is integer, then 

(-1) 2 / 3 = (-IY W = (_l)2/Hl/2-aarg(z)/(27r)J _ ^ 

so L can be simplified to 1. 
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Remark 12. In the next most common case where /3 is a half-integer, L\ can be expressed 
more candidly as 

_ 1 if mods(7r-aarg(z), 2vr) > 0, 



-1 otherwise, 



and L 2 can be simplified to 



L 2 = (_i)W(2-) J X if ~ n < ar s( c (*)) + arg(2 a ) < 7T Q 

—1 otherwise. 

1 if — 7r < arg(c(z)) + arg(z Q ) + arg I 1 + — ^ — I < n 



c(z)z a J ~ ' (70) 

1 otherwise. 

It is easy and worthwhile to test for these cases. Alternative (169 j) is preferable when w is 
a constant. 

Techniques described in subsection 12.31 can further simplify local equivalents to L 2 
and Z/2 _ particularly for real z. 

Also, if the local equivalent simplifies to a piecewise constant that is one value £q at 
z = for which 0£o simplifies to 0, but another value C* everywhere else, then we can 
simply use £* for all of the terms having positive degree, because those powers of z are at 
z = 0. This can dramatically simplify the result because often there are no non-positive 
degree terms or only one. For example, with real x there is only one non-positive degree 
term for 



series 



((-1 + x -ix 2 ) 7/2 , x = 0, o(x 2 ) ) -)■ < ' j ix-\ ix 2 + o(x 2 ). 

V ' ' ' V ') \\i otherwise^/ 2 8 K ' 



Otherwise if L doesn't simplify to a constant, then we have to include a perhaps 
complicated factor multiplying every term of the result, because extracting the dominant 
coefficient and exponent are mandatory for numeric powers of series. This is significantly 
more annoying than the possible one piecewise-constant term for a logarithm of a series. 
If we choose to distribute non-constant L over every term, then the result is significantly 
bulkier and less intelligible. If instead we choose to return an undistributed product 
of this factor with a sum of terms, then we have denied the user what they implicitly 
requested by using a function named series - a sum of terms. 

If we are implementing the full generality of a hierarchical series, then such recursively- 
represented series have some appeal, because the jumps in the rotation angle make it have 
some properties of some essential singularities, which dominate any power of z. However, 
one can argue that such jumps belong in the coefficients because the magnitude of the 
rotation factor is always 1, whereas the logarithmic singularities that we already allow 
in the coefficients have infinite magnitude at z — 0. 

It is easy for a human user to use an expand (...) function to distribute the factor over 
the terms after seeing an unexpanded result. However, the user might be other functions 
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whose authors must be knowledgeable enough to realize that they should always apply 
expand (. . . , z) to the result of series (...) because it could be an expandable product. 

If we are not implementing the full generality of hierarchical series, then a more serious 
disadvantage of not automatically distributing the factor is that it requires a special 
field for a multiplicative factor in the series data structure; and one such specialized 
exceptional field is inadequate for adding series containing such multiplicative factors if 
the series have different dominant exponents. 

Expressions such as (— l) 2 ^ can alternatively be expressed as L := e 2 ^ 7 ™ or repre- 
sented that way internally. However for presentation in results, (— 1) 2/3 avoids a perhaps- 
unnecessary % and more obviously indicates that the factor has magnitude 1. 

Remark 13. Whenever z is real and L contains more than one piecewise factor, their 
product can often be combined into a single factor by separately simplifying their prod- 
uct under the alternative constraints z < 0, z = 0, and z > 0, then forming a single 
piecewise or constant factor accordingly if the three values are constants. This process 
often simplifies the product to the constant 1 when (3 is a half integer. 

3.2 Real-branch series of fractional powers 

For computing a fractional power of a series we must fully distribute the exponent over 
the dominant coefficient and dominant power, then combine the two exponents of the 
latter into a product. A way to accomplish this correctly for real-branch mode is to use 
in this mode the following rewrite rules for reduced exponents for a real expression t 
together with integers m and n: 

arg ( t 2m/(2n + l)) ^ Q, (71) 

arg ( t (2m+l)/(2n+l)) _> arg (t) . (73) 

4 Branch bugs for inverse trig & hyperbolic series 

Kahan [3] is often cited as a standard for defining the branch cuts of the inverse trigono- 
metric and inverse hyperbolic functions, together with the principal values on those 
branch cuts and on the branch points that end them. In comparison to Derive, his 
definition of arctan (...) has the disadvantage of violating the useful identity 

arctan I — ^=^= J = arcsin (z) . (73) 



vT^ 



However, this section uses Kahan's definitions, which are used by TI-Nspire. 

This section doesn't discuss the six secondary functions such as arcsec (z) := arccos (1/z) 
and arccot (z) := n/2 — arctan (z) because their series are easily computed from such 
definitions. Beware, however, that different mathematical software might use different 
definitions. For example, Mathematica uses 



vr/2 if z = 0, 

arccot (z) := ■{ % 

— (In (1 — i/z) — In (1 + i/z)) otherwise, 
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which is not equivalent everywhere to ir/2 — arctan (z). 

One way to compute series for the primary inverse trigonometric and inverse hyper- 
bolic functions is by their integral definitions. However, the handling of non-constant 
coefficients that can include nested logarithms and piecewise constants is then problem- 
atic, and this method alone doesn't specify the degree-0 term of the result. 

Consequently, this article instead uses Kahan's definitions in terms of logarithms and 
powers. 

Many of the tests in this section entail comparing the imaginary part of a series V with 
0. For both real and complex z these tests can often be simplified by determining critical 
angles and perhaps also computing I C (V,8 C ), as described in Sub-subsections 12.1.11 and 
12.1.21 Some of the following formulas instead entail comparing the real part of a series V 
with 0. They often can similarly be simplified by determining critical angles and perhaps 
also computing I c (iV, 9 C ), which maps the real part to an imaginary part. 

4.1 Series for arctanh 

The inverse hyperbolic tangent of a series U can be computed from the identity 

ky(l + U)-\n(l-U) 
arctanh (U) = — -. (74) 

The logarithms in formula ( I74"|) can contribute expressions involving u, (f) or ip as 
discussed in Section [21 The result is candid in most cases. However, if series U = 
cz a + g(z) has a negative dominant exponent a, then the zero-degree term computed by 

flUD is 

In (cz a ) - In {-cz a ) + (ft (U + 1) + ft (-U + 1)) . (75) 

This expression can and should be replaced with the simpler local equivalent 

{±1 if = 0, 
1 if 9f (U) > V Of (U) = A SR (U) < , (76) 

— 1 otherwise, 

which is also less prone to rounding errors. 

4.2 Series for arctan 

The inverse tangent of a series U can be computed from the identity 

arctan (U) = —i arctanh (ill) . (77) 
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The degree-0 term of arctan (U) is 

if a > 0, 

±tt/2 if a < A z = 0, 

tt/2 ifa<0A(R(u)>0Vfi(u) = 0AS (u) > 0) , 

-tt/2 if a < A (»(m)<0V»(u) = 0AS(h)<0), 

< | (-« In (^) + Q [ibz a + i/t(z))) if U = i + bz° + h(z), 
|(zln(=^) -Sl(-ibz?-ih(z))) iff/ = -i + bz' 7 + h{z), 
i arctanh(t) + |^ >0 (1 - t + ig{z)) if U = ti + g(z) A t > 1, 
i arctanh(t) - -fl^o (1 + t - ig(z)) if U = ti + g(z) A t < — 1, 
ar ct an ( c ( z ) ) ot herwise . 

(78) 
Formula (ITT)) should automatically achieve these simplified forms except probably for the 
three cases where a < 0. 

Table [3] lists some simplified correct 0-degree terms for arctan (. . .). Corresponding 
examples for arctanh (. . .) can be derived from f ITTj) . 

4.3 Series for arcsinh, arcsin and arccos 

The inverse arcsine, arccosine and hyperbolic arcsine of a series U can be computed from 

arcsinh (U) = \n(u + y/l + uA , (79) 

arcsin (w) = —i arcsinh (iw) , (80) 

7T 

arccos (w) = — — arcsin (w) . (81) 

The square root might contribute a factor containing expressions involving arg (...) to 
all of its result terms. The subsequent logarithm in identity (1T91) might then contribute 
horrid nested instances of arg (...) to a 0-degree term in the result. 

Mathematica 7.0.1.0 avoids this by instead using correction terms and factors that 
are specific to inverse trigonometric and inverse hyperbolic functions: The resulting ar- 
guments of arg (. . .) are unnested and entail the terms of U, which are often fewer and 
simpler than those of 1 + U 2 and U + \[\ + U 2 . Extensive numeric experiments provide 
convincing evidence that they are correct in most instances. Most of the formulas in this 
subsection and the following one are simplified and corrected versions of those correc- 
tions. When using these corrections and identity (jT9j) . it is important to use the rewrite 
In (pz a ) — y ln(p) + aln(z), but suppress the branch corrections for logarithms and frac- 
tional powers discussed in previous sections. It is also important to compute intermediate 
series such as a/1 + U 2 to sufficient order so that the final order is as requested. 

A correction term or factor is unnecessary for arcsinh (U(z)) if the dominant degree 
of U(z) is positive or if the dominant degree is and U (0) isn't on either branch cut or 
either branch point. Otherwise, here are the different cases: 
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4.3.1 Case arcsinh (U(z) = ci + g(z)) with c < — 1 V c > 1 

If the dominant term of U(z) is ci with c < — 1 V c > 1, then 

msignfc) i / n\\ 

arcsinh (a + <? (2)) = \- T ■ (arccosh (c) + o [z J J 

where 

1 if sign(c)3*(^)) > 0, 
•1 otherwise. 

4.3.2 Case arcsinh (U(z) = ci + g(z)) with c = 1 V c = — 1 



(82) 



If the dominant term of U(z) is ci with c=lVc = — 1 and the next non-zero term is 
bz a , then 



where 



arcsinh (ci+te'+Zi^)) = — + (-l) P Q ■ (icV2ibz a/2 + o (W 2 ) 

ich(z 



I arg^c6+ ^ o 



o" arg(zj 



g 



ich(z 



2n 



■1 ifarg(6) = ^A»^)<0, 

1 otherwise. 



Expression icft H )■ ic6 as z — > 0, so the critical angles for P are 

(2ra + 1) 7r - arg (icb) 

— IT < U r — < 7T 



<7 



for all integer n satisfying 

arg {icb) 1 a 



2n 



2 2 



< n < 



arg (ic6) 1 a 

2^ 2 + 2 



There are no such angles if < a < 1 and arg (icb) is sufficiently close to 0: If 

(a — 1) n < arg (icb) < (1 — a) n, 

then P = 0, which is easily tested. 

If < o" < 2, then a more candid and accurate representation of (— l) p is 



I if — n < arg ( icb + 
-1 otherwise. 



ich(z) 



+ aarg(^) < 7r, 



(83) 



M) 



J5) 



(86) 



57) 
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4.3.3 Case arcsinh (. . .) with negative dominant exponent and imaginary dom- 
inant co efficient 



Use formula (|T9|) including the angle rotation factor for the square root and the correction 
term for the logarithm. Brace yourself for a truly ugly result. 

4.3.4 Case arcsinh (. . .) with negative dominant exponent and non-imaginary 
dominant coefficient 

If the dominant exponent a of U(z) is negative and the dominant coefficient c isn't pure 
imaginary, then 

arcsinh (cz a + g{z)) = (-1) w N ■ UitW + ln ^ ^ + a In (z) + o (z°) ) . (89) 



where 



W 



N 



arg c + 
1 \ V 



+ 2a arg(z) 



2n 



1 if K(c) > 0, 
■1 otherwise. 



(90) 



(91) 



Expression I c H ) — > c 2 as z — > 0, so the critical angles for W are 



(2n + l)7r-arg(c 2 ) 

7T < &r = : < 7T 



2a 



for all integer n satisfying 

arg (c 2 ) 1 a 
2tt 2 ~ 2 



< n < 



arg (c ) 1 a 



2tt 2 2 



(92) 



(93) 



There are no such angles if —1/2 < a < and arg (c 2 ) is sufficiently close to 0: If 

(a-l)7r < arg(c 2 ) < (l — a)ir, (94) 

then W = 0, which is easily tested. 
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If — 1 < a < 0, then more candidly and accurately, 

2 N 



w 



1 ifarg|(c + ^-) j + 2a&rg(z) > vr, 
-1 ifargj(c + ^) ) + 2a arg(Y) < -vr, 



otherwise; 

'-1) W = { -o.y-. z a J J • °W-" ^ 




-1 otherwise. 



Some incorrect series for arcsinh (...), arcsin (...) and arccos (...) in current imple- 
mentations are attributable to incorrect or omitted correction terms and/or factors. Ta- 
ble H] lists some examples for arcsinh (...). Corresponding examples for arcsin (...) and 
arccos (...) can be derived from flHUl) and (1811). 



4.4 Branch corrections for arccosh 

Subsection 14.31 expressed the series for arcsin and arccos but not arccosh in terms of the 
series for arcsinh. The inverse hyperbolic arccosine of a series U can instead be computed 
from either of 




arccosh([7) = 2 In I yj —^- + yj —^- I , (97) 

arccosh = ln((J + y/U - ly/U + l) . (98) 

The latter is slightly slower asymptotically because of the series multiplication. The 
alternative that gives a simpler result might depend on the general expression simplifier 
and the case, such as whether or not the dominant degree of U is negative. 

Each square root in these formulas might contribute a factor containing expressions 
involving arg (...) to all of its result terms. The subsequent logarithm in identity (1971) or 
( 198]) might then contribute horrid nested instances of arg (...) to a 0-degree term in the 
result. Here are some simplified correction terms and factors that avoid this: 

Let the dominant term of U be cz a , and let the sum of the remaining terms be g(z). 
Let bz a be the dominant term of g(z), and let h(z) be the remaining terms of g(z). A 
correction term and/or factor is unnecessary if a = A (S (c) 7^ V c > 1). Otherwise 
we have the following cases: 

4.4.1 Case arccosh (U( z)) with dominant deg & dominant coefficient < -1 

If a = and c < — 1, then 

arccosh (c(z) + g (z)) = 2in J + arccosh (c(z)) + o (2 ) (99) 
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where 



J 



-1 if9(<7(*))<0, 

otherwise. 



(100) 



4.4.2 Case arccosh (U(z)) with dominant degree and dominant coefRcient in 
(-1, 0) or (0,1) 



If a = A (-K c < V 0<c<l) then 

arccosh (c + g (z)) = K ■ (arccosh (c) + o (z )) 

where 

, (-1 if%(g(z))<0, 



k = (-iy 



1 otherwise. 



:ior 



(102) 



4.4.3 Case arccosh (U(z)) with dominant degree and dominant coefRcient 1 

If instead a = and c—1, then 

arccosh (I + bz a + h (z) ) = 5£ • (V2bz a/2 + o (W 2 ) V (103) 



where 



B 

E 

h{z) 



arg [b + 

1 V z 



Hz) 



+ a arg (z) 



2n 



-1 if axg(6) = 7r A A ( ^) <(). 
1 otherwise. 



Expression b -\ — ^ — )• b as z — >• 0, so the critical angles for B are 

(2n+l)vr-arg(6) 

— 7T < tA = < TV 



< n < 



arg (6) 1 a 
~1^ 2 2 



for all integer n satisfying 

arg (b) 1 a 
_~1>^ 2 ~ 2 

There are no such angles if < er < 1 and arg (b) is sufficiently close to 0: If 

(<7 — 1) 7r < arg (6) < (1 — <r) 7T, 

then i? = 0, which is easily tested. 

If < a < 2, then a more candid and accurate representation of B is 

1 a arg (2) < n, 



1 if — 7T < arg ( b + 
-1 otherwise. 



(104) 
(105) 



(106) 



(107) 



(108) 



(109) 
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4.4.4 Case arccosh (U{z)) with dominant degree and dominant coefficient -1 

If a = and c = —1, then 

arccosh (-1 + bz a + h (z) ) = iirC + CBE ■ (iV2bz a/2 + o (z a/2 ) 



where 



C 



1 if %(bz a + h(z)) < 0, 
-1 otherwise. 



110) 



4.4.5 Case arccosh (U(z)) with negative dominant degree 

If instead a < then 



arccosh (c,2 a + # (z)) = 2m (D 1 + M) + In (2c) + a In (2) + o (2 ) (111) 



where 



£>i 



A/ 



1 ar n c+ ~ 



a arg (z) 



2^ 



1 if arg(c)=7r AS(^|<0, 
otherwise. 



9{z) 



'1121 



113) 



Expression c + ^- — > c as z — > 0, so the critical angles for £>! are 

(2n + l)7r-arg(c) 

— 7T < C7 r = < 7T 



O 



for all integer n satisfying 

are; (c) 1 a 



2tt 2 2 



< n < 



arg (c) 1 a 
~~ 2^ 2 + 2 



There are no such angles if — 1 < a < and arg (c) is sufficiently close to 0: If 

— (1 + a) n < arg (c) < (l + a)7r, 

then D\ = 0, which is easily tested. 

If —2 < a < 0, then a more candid and accurate representation of D\ is 

( g( z )\ 

1 if — 7r < arg I c H — J + a arg (z) < n, 

— 1 otherwise. 



(114) 



115) 



(116) 



117) 
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4.4.6 Case arccosh (U(z)) with positive dominant degree 

If instead a > 0, then 



arccosh (cz a + g (z)) 



■i) D n-v D2 G-{^+o(z°) 



in 



'1181 



where 



Do 



1 1 

2 ~ 2^ 



/ 



/ 



arg 



V 



V 



c + 



9{z) 



a arg(zj 



(119) 



G 



-1 if z ^ A 3 (c) < V 

arefc) = tt A $ ( ^ I < V 



are; (c) = A S 



»(*) 



cz a + #(z) 



(120) 



<0, 



1 otherwise. 
Expression —1/ (c + g(z)/z a ) — > — 1/c as z — >■ 0, so the critical angles for Z} 2 are 



a 



for all integer n satisfying 

arg (—1/c) 1 a 



< n < 



arg (-1/c) _ 1 a 
2tt 2 2 



2tt 2 2, 

There are no such angles if < a < 1 and arg (c) is sufficiently close to 0: If 

(a — 1) 7r < arg (—1/c) < (1 — a) tt, 

then D 2 = 0, which is easily tested. 

If < a < 2, then a more candid and accurate representation of D 2 is 



(122) 



(123) 



1 if — 7r < arg ( —1/ ( c 
— 1 otherwise. 



9W 



+ a arg (z) < n, 



(124) 



When using these corrections and identity (197|) or (198]) . it is important to use the 
rewrite In (cz a ) —¥ ln(c) + aln(z) but suppress in these formulas the logarithmic and 
fractional power adjustments discussed in Sections |2] and [3j Also, sub-expressions such 
as 

-i - g(z) 



c + 



9{z) 



and 



cz a + g(z) 



29 



in (I119p and ( |120p can and should be approximated by appropriate-order series. Note 
that a truncated series for these two sub-expressions usually can't be exact for non-zero 
g(z). Therefore we should expect some narrow cusps in which the corrections are not 
piecewise constant. If g(z) is known to o (z m ), then we can compute the series for these 
two sub-expressions to o(z m ~~ a ). 

For both arcsinh (U(z)) and arccosh (U(z)), if we are unable to narrow the choice to 
one of the above cases for reasons such as c containing insufficiently constrained indeter- 
minates other than z, then we should construct a piecewise result from all of the cases 
that we can't preclude. 

Table [5] contains some correct results for arccosh (. . .) series. 

Summary 

The generalization from Taylor series to generalized Puiseux series introduces a surpris- 
ing number of difficulties that haven't been fully addressed in previous literature and 
implementations. The most serious of these is incorrect results for expansion points that 
are on a branch cut or a branch point. Formulas are presented here that correct this 
for logarithms, fractional powers, inverse trigonometric functions, and inverse hyperbolic 
functions. These corrections typically entail an additive piecewise-constant multiple of 
27T2 and/or a unit-magnitude piecewise-constant multiplicative factor. Some of these cor- 
rections are applicable even to Taylor series. There are many alternative formulas for 
these corrections. The alternatives presented here are chosen: 

• to reduce catastrophic cancellation when evaluated with approximate arithmetic, 

• to candidly reveal the piecewise constancy, 

• to reveal the boundaries of the pieces as explicitly as is practical, 

• to simplify significantly where practical, such as for real expansion variables or 
numeric coefficients. 

• to be as concise as possible subject to the above goals. 

Tests, formulas and algorithms are given that compute simplified versions of these cor- 
rection terms and factors near the expansion point for both real and complex expansion 
variables. 
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Table 1: O-degree term of series (lnw, var = 0, o (var n )) with n>4, x, y G K, z = x+iy: 



Alternative O-degree terms of In u near variable = 



If 



u 



why 



la 



z 2 + z 3 



2m if 3( 2 )<0A5R(z) > %(zy/2 + --- + o(<S{z) 
ln(z 2 ) + { -2m if 3(z)>0A5J( 2 ) > 3 (z) 2 /2 + ■ ■ ■ + o (3 (z) 







otherwise 



(TO 

(E3D 

rem. [6] 



lb 



2\nz + 2m 



7i — 2 arg z 

2m 



2ztt, 3(z)<0AK(z)>3 (z)72 + 
-2m, 3(z)>0A3J(2)>3 (z) 2 /2 + 
0, otherwise 



(USD 

rem. [6] 



2a 



z 2 + z 3 e 2 



2i7r if &rg(l+z + z 2 + ho (z n )) + arg(,2 2 ) < -n 

ln(z 2 ) + I -2iir if aig(l+ z + z 2 ^ ho (z n )) + arg(z 2 ) > vr 







otherwise 



2b 



2\nz+27ri 



7r — 2 arg z 



2jt 



2m, arg(H h o(z n )) + arg(z 2 ) < — n 

+ ^ -2m, arg(H h o(z n )) + arg(z 2 ) > n 

0, otherwise 



2c 



lnfz 2 ) + (arg (z 2 + ho (z n )) — arg(l+ ho (z n )) — arg(z 2 )) i 



2d 



21n(z) + (arg(> 2 H Vo(z n )) - arg(H \-o(z n )) -2arg(z))i 



. Z -W- Z V* 



ln(- 



-7/6 



) or —7 In (z) /6 



ITT 



if 3 (z) > 



— 7tz otherwise 



if (S(z) >0AK(z) <3(z) 2 /2 + - 
(9f(z) >0AS(z)> 3?(z) 2 /2 + -' 
— i7r otherwise 



prop [5] 
or (ED 



27T 



-1-Z 2 -Z 3 



• + o(3(z) n ))v 



(EHD 

rem. [6] 



-l-z 2 e 2 



m if$(z 2 + 2z 3 + z 4 + --- + o(z n ))<0 
-in otherwise 



-l + zz 1/4 



/7T 



■1—izi + ^ 



27T if = 

-in otherwise 



-1+iz* +z 



m if z < 
— «7r otherwise 



rem. [8] 



C + 2T 



In (z 2 ) if c = 

In (c) - 2in if arg (c) = n A S (z 2 ) < 

In (c) otherwise 



rem. [9] 



10 



-x + e 



In (x 



prop[3] 



11 



x 2 + x 3 e x 



In (x 2 ) , or 21n(|x|), or 21n(x) + 



-2in if x < 
otherwise 



x 4/3 _|_ x 2 



In fx 4/3 



prop [3] 
or flolj) 



12a 



rem. [7] 



12b 



real 
branch 



ln(x 4/3 ), or 41n(|x|)/3, or 41n(x)/3- 



4ivr/3 if x < 







otherwise 



sec. 
12X21 
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Table 2: For n>4, x,yE.M., z = x+iy: series(u /3 , var = 0, o(var n )) —¥ (—1) A (£?+•• 



//■ 



u h 



A, or L = (-1 



5 



why 



la 



{z 2 + z 3 f 2 



L-- 



1. (3( 2 )>0AK(z)>3(z)72 + . 
V (Q(z)<0A &(z) >3(z) 2 /2 + . 
— 1, otherwise 



-o(3(*) n )) 
-o(9f(z) B )) 



(HQD 

rem. [6] 



lb 



I if (3£(z) > V £(z) = A 3(z) > 0) 
= (— 7r < arg (z 2 ) + arg (1 + z) < 71") 
•1 otherwise 



(EHD 

cam 



lc 



A = ^ (arg (z 2 + z 3 ) — arg (1 + z) — 2 arg (2;)) 



2a 



(z 2 — iz 3 ) 



3N 3 /2 



L 



1 if ft (2) > V 
-1 otherwise 



A 3 (2) > 



dSQD 



2b 



A= LV2-(axgz)/2j 



3a 



fz 2 + ^ 3 ) 7/4 



L=< 



X (3W>0A^)>3(2)72+...+o(3W)) 

v(3(z)<0AK(2)>3(2)72+...+o(S( Z ) n )) 

(3(z)>0Afi(2)<3(z)72+...+o(3(2)")) 



y 7/2 



-i, otherwise 



^72" 



(EZD 

rem. [6] 



3b 



A 



(- 1+lz y^ +z f / 



An 



(arg (z + z ) — arg (1 + z) — 2 arg (z)) 



L= 1 



—% 



W 



(J3T1I52D 



(-l+izV 2 +z) 



L= 1 



rem. [8] 



6a 



' 1 2 3\ 3 / 2 

— 1 — Z — Z ) 



-1 if (y > A x < y 2 /2 - 3y 4 /8 + 
L = ( (y < A x > y 2 /2 - 3y A /8 + 

1 otherwise 



-)V 



—% 



dSHD 
dZQD 

rem. [6] 



6b 



A = 2^ ( ar ^ 



-l-z"-z d 



arg^ + z^ + z' 1 ) — 7r) 



(c+^ 2 ) 3/4 1 c^O 



/ if arg c = 7r A S (z 2 /c) > 
1 otherwise 



.3/4 



rem. [9] 

am 



(-z-^) 1/2 



L 



-1 if arg (z) < 

I otherwise 



iz 



(ln(z) + z) 3/2 



rem. [2] 

am 



L = 1 



111(2) 



10 



[iz + z 



2\ 3 h 



L 



-1 ifK( 2 )<0A3(z)>0 
I otherwise 



(-1) 4 zi 



(1231) 
(El 



11 



>- : + 1) 



3/2 



•1 if z < 

I otherwise 



3/2 



(-2 + a;) 3/2 



prop. [5] 
(El 



12 



L 



prop. [3] 



13 



(-1 " Vx) ' 



/2 



1 if x > 
•1 otherwise 



14a 



(x 4 / 3 + x 2 ) 



3/2 



I 



■1 ifx<0 

I otherwise 



xr 



( I23T691 



14b 



real branch 



L 



;r 
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rem. H] 



Table 3: O-degree term of series (arctan (u) , z — 0, o {z n )) with zGC, iGK, n > 3 : 



# 



» 



A correct O-degree term near z = 



2T 2 + 2T 1 



7T 



±1 if 3 = 

1 if - tt/2 < arg (^~ 2 + 2" 1 ) < vr/2 
-1 otherwise 



I »*(* + *' + £ + 



2z + ze 2 



In 3 



+ o U r ' 



<0 



7T 



otherwise 



~tc In 3 

- H z 

_2 2. 



2z + ^V 



In 3. I 2 if ^° 

2 I otherwise 

2 



-1i + ze z 



~7T In 3 



-% - z l ^e z 



i + ze 2 



7 ? ? 

vr - In (zz) + In 2 10 if arg [l+z-\ 1 1 h of/ 

-z — ^ \ 2 6 

7r otherwise 



arguz) < n 



i + ize 2 



ln(- 
2. 



1 

_2_ 



zln 



— z + z 2 e 2 



— iz 



( / ^2 T3 

7r if arg I 1+2 + — + — + 
2 6 



+ < 



vr if arg( 1+2 + — + — - 
2 6 

otherwise 



+ o(z n ) I + arj 
+ o(z n ) ] + arj 



— ZZ ) < — 7T 



"ZZ > 71 



i In (z/4) 



— z + iz^ 2 e z 



10 



x 2 + x * 



7T 

2 



±1 ifx = 

1 if x > 
-1 otherwise 



- if x > 

2 I otherwise 

2_ 



11 



2z + x 3/ V 



In 3 



12 



2i + W r t 



7r/2 + z(ln3)/2 



- if x > 

2 I otherwise 

2_ 



13 



-2z + a; 3/ V 



In 3 



14 



— 2i — x 3/4 e :c 



— 7r/2 — i (In 3) /2 



- In (ix 5 / 4 /2 



15 



z + x 5 / 4 e x 



16 



-z + x 5//4 e a: 



In (-ix^ 4 /2) 
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Table 4: First two non-zero terms of series (arcsinh (it) , z = 0, o {z°°)) with z G C, x G 



u 



First two non-0 terms of arcsinh u near variable = 



Sec. 



-2% + z 2 + z 



2 i „3 



27T 



1 if » ^ + z 3 > \ / , ~ z',2 2 

i l , In (2 + V3) + -= + 

•1 otherwise / V V3 



axu 



2z + Z 2 + ,2 



2 _i_ „3 



Z7T 

T 



1 if^(z 2 + ^ 3 ) >0W In (2 + ^3) iz 2 



-1 otherwise 



v/3 



iXU 



T74 

y + ln(2 + v^)-^ + . 



2i + ^ 1/4 



IHO 



z + z 2 + z 3 



T + M > 



1 arg (—1 + iz) + 2 arg 2: 
2~ 



2n 



[iz + 



l4~3T2l 



i + iz 2 + iz 3 



- + (-1) 

2 v y 



1 arg (— l+iz)-|-2argz 
2" 



2^ 



1" arg (—1 + iz) + 2 arg z 

2~ 2^ 



1 ^(z)>0 \ 
■ 1 otherwise J 



2z+---) 



BH 



-i + iz 2 + zz 3 



Z7T 



(iV2z+---) 



I4T3T21 



-i — iz 2 + z 3 



-jH-D 



T arg (— l-Mz) + 2arg,z 
2" 2^ 



(fl, R(,)<0 \ 
V I — 1 otherwise I 



,+ •■■) 



IQ31 



iz" 1 + 1 



ugh! 



14331 



1" arg(l + 2z + 2 2 )+4ar^ 
2 



9 

10 

11 

12 

13 

14 

15 

16 

17 

18 



z~ 2 + z~ l 



-l) w (\n2-2\nz+7iiw + z + - ■ 



w ■ 



2n 



B33 



-2% + x 3/2 



7T7 



1 if x < 



— -ln(2 + v / 3) . 
2 1—1 otherwise 



z.r 



3/2 



v/3 



+ 



iXU 



372"" 



-2% - x 3/2 



7T2 



, ■n(2 + V3)j+- 7r + 



2.7; 



Hxn 



2z + x 3/4 



7T? 

y 



, , ^ 1 if x > \ ix 3 l A 

+ In 2 + v^) < " pr- 

v 1-1 otherwise/ v 7 ^ 



+ ■■ 



I4XT1 



T + ^ + ^-JJ 



2z + x 3 / 2 



I4XT1 



z + x 



3/2 



7TZ 



if x> 

•1 otherwise 



((1-1)^3/4 + ...) 



l4~3T2l 



i — x 3//2 



7T/ 



__(l + ,) x 3/4 + ... 



IHOl 



-i + x 



5/2 



7TZ 



if x > 

-1 otherwise 



((l + z ) x 5/4 + ...) 



I4~3T21 



-i + x 



3/2 



7T/ 



^ + (1 + z) x 3 / 4 + 

w 



l4~3T2l 



9 1 

a; + x 



,0 if x > \ 

In 2 -21nx + { ~ +x + 

2iri otherwise / 



MM 
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Table 5: 1 st two non-0 terms of series (arccosh (u) , z = 0, o (z°°)) with z G C, x G 



1 
T 

3 



u 



1 two non-0 terms of arccosh u near z = x = 



Sec. 



-2 + z 2 + ^ 



2 i ^3 



In (2 + V3) 



/" if3(z 2 + z 3 )>0 
-Z7r otherwise / \/3 



+ 



EH 



-2 + iz 1 / 4 + z 



(In (2 + yg) + in) - iz 1 / 4 /V3 + - 



su 



- + z 2 + z 3 



I if 3 (z 2 + z 3 ) > 
-1 otherwise 




I4A21 



1 + z 2 + z 3 



1 if — 7r < arg (1 -+- z) + 2 arg z < 7r 
-1 otherwise 



^ + 7! + 



I4A31 



■1 + z z + z 



2 i ^3 



1 3 (z 2 + z 3 ) > 
-1 otherwise 

\B = 



(in + B- {-i^2z + ■■■)) 
1 arg (1 + z) + 2 arg z 



1) 



2tt 



I4A41 



z~ 2 + z" 1 



2ivr 



7r — arg (1 + z) + 2 arg (z) 

2^ 



+ In 2-2 In z + z + ■ 



MB 



■1 otherwise / V 2 



z + z^ 



I A 

AD 2 = 



arg (1 + z) + argz 

2^ 
arg (— 1 + z + z 2 + ■ 



argz 



2n 



MM 



-2 - x 3 / 4 



In (2 + V3) 



,3/4 



27T 



4 if x < V ^ 
otherwise / y/3 



EO 



_9_ 
10 

IT 

12 
13 

14 



-2 + x 



(In (2 + yg) + in) - ix/V3 + ■ ■ 



SU 



--VX 



I if x > 
-1 otherwise 




2JX 1 / 2 
73 



I4A21 



X /2 + X 



ni 



/3 - 2ix/V3 



Q31 



1 — x — x 3 



I if x > 
-1 otherwise 



\/2ix -\ =x 2 + 

V2 



I4A31 



1 — x 2 — x 3 / 2 



v / 2ix + ix 3 / 2 /\/2 + 



Em 



•l + x 2 



7ri + 




if x > 
otherwise 



(-v / 2^ + ---) 



MM 



15 



— 2 — ] 



-2 In x + In 2 



ni if x > 

3ni otherwise 



3ix 



+ 



-TJT 



2x — x 3/2 



-2 In x + In 2 + ni) - x 3/2 H 

7ri/2 - 2ix H 



Q31 



16 

17 



-x 2 + x 



14X51 
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